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Let an elastic plate occupy a region D with boundary C. On a portion 2 of C the 
deflection, slope, moment, and shear are measured. In terms of these data, and crude 
bounds for the maximum deflection and moment, pointwise bounds for the value of the 
deflection and its derivatives in D are obtained. These bounds are such that if the 
data on 2 and the loading perpendicular to the plane of the plate tend to zero, the deflec- 
tion tends to zero at every point in D. Similar bounds are obtained in the slow flow 
problem for a viscous fluid. 

1. Introduction 

The following plate problem is considered. Let an elastic plate occupying a region D 
in the xy-pb&ne be subjected to a load p(x } y) normal to the plane of the plate. The deflection 

■v satisfies the equation 

A2 /d 4 v , _ d 4 y . d 4 A irj> n ,. 

A H^ +2 ^w + ^ 4 h p/ ^ (L1) 

where D denotes the plate rigidity (see [3, p. 488]).' We assume that on a portion 2 of the 
boundary of D we are able to measure the deflection, slope, the sum of the two principle mo- 
ments, and the normal component of shear with known error; or cquivalently we assume 



f (v-f¥d<r<ai, f («. t -F t )(v, ,-F t )d<r<a 2 

f (Av-gyd<r<a 3 , f (A» { -G,)(Av, t-G f )d<r< ai 



(1.2) 



where/, F i} g, and G f are the measured values of v,v ii} Av, and Av, ( respectively. In (1.2) 
the comma denotes differentiation and a repeated subscript i implies summation over i (i 
= 1, 2). The A symbol is used throughout to denote the Laplace operator. 

In (1.2) then the a t are assumed to be small known constants. We wish to obtain upper 
and lower bounds for the solution v at any point in E. 

It is well known that as it stands the Cauchy problem for the biharmonic equation is 
not well posed (see Hadamard [1]). The solution does not depend in a continuous way on 
the boundary data. We leave aside the question of existence of solution and merely remark 
that in most physical situations in which one would be interested in considering such a prob- 
lem the question does not arise. For instance, we might wish to treat an ordinary plate 
boundary value problem, but find that part of the boundary is inaccessible for measurement 
•of boundary data whereas on another portion of the boundary it is possible 4 to measure v, 
bv, Av, and _d^ (Av) . 

If in addition to the Cauchy data we have uniform bounds for both v and Ay throughout 
D then the solution to the plate problem will be seen to be stable. This is in fact an im- 
mediate consequence of the results of Laurentiev [2]. 
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In order to simplify the procedure somewhat we decompose the solution v into 

v=v p +v', (1.3) 

where v p is any particular solution of (1.1) defined in D and v f satisfies 

AV = (1.4) 

in D. We approximate v' by a biharmonic function <p and set w—v' — <p. Thus w satisfies 
u\D 

A 2 w=0. (1.5) 

Also bounds for the Cauchy expressions 



rji= W 2 da, rj 2 = W ti W ti d(T 

6i= (Aw) 2 d<r, e 2 = Aw,iAw ti d<r 



(1.6) 



are known. In fact 



Vi~- 



- f i(v'+v P -f)-(<t>+v P -f)Yda<2a 1 + 2 f {<t>+v P -j)H° (1-7) 



with similar inequalities for rj 2 , e h and e £ . It is clear that if a\ is small and can be chosen to 
approximate/-^ in mean square on 2 then ^ will be small. A similar statement holds for rj 2l 
e h and e 2 . We shall show in this paper that it is possible to obtain an a priori bound for 
w of the following type at any point in R: 

|^(P)| 2 <^ 1 M!-^{a 1 r ?1 +a 2 r 72 }^+yl 2 ^- 5 n 6iei + 6 2€2 } 5 2 7 ( L8) 

where A h A 2 , d\, a-:, b h b 2 , are determined constants, <5i and 8 2 are constants satisfying* 

0<8,<8i<l, (1.9) 

and Mi and M 2 are the uniform bounds for w, and Aw in D. It is clear that knowledge of uniform 
bounds for v' and Av' implies knowledge of uniform bounds for w and Aw since w-v' is a known 
function. 

With an inequality of type (1.8) it is clear from (1.7) and the corresponding inequalities 
for rj 2 , €i, and e 2 that if the a t are sufficiently small and the Rayleigh-Ritz technique is used to 
make the integral expression on the extreme right small then we obtain close upper and lower 
bounds for v'(P). This paper will be concerned then with the establishing of (1.8). 

2. A priori Bounds 

We consider a two-dimensional region D bounded by a closed surface C. On a portion 
2 of C, Cauchy data are prescribed. We_denote the remainder of the surface (C-2) by 2. 
(As will become apparent later the portion 2 need not actually be prescribed.) We now choose 
the origin of a polar coordinate system at a point P which has the property that there exists a 
circle K(r) of radius r and center at P such that the intersection of C and Kir) contains no 
point of 2. (If two or more disjoint regions are formed by the intersection of Kir) and D, it is 
sufficient that the boundary of any one of these regions contains no points of 2.) We denote 
by D r the intersection of K{r) and D (the portion whose boundary contains no points of 2 if 
two or more disjoint regions are formed) and require that P lie outside of D r . We denote by 
R the radius of the largest circle satisfying the above conditions (for a fixed P ) and by r the 
radius of the smallest such circle. As indicated in the previous paper [5] we locate P in such a 
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way that K(r ) is tangent to D at some point of 2. Finally we denote by co(r) the portion of 
the surface of K(r) which lies interior to C. 

As in [5] we introduce for any function u the notation 



u 2 do 



where da denotes an element of surface area on 2, and let 



(2.1) 



^>=U[(§)'+(!)>* 

Dr 



In [5] we proved that if u satisfies 



then the functional 



Au=0, \u\<MinD, 



v(r 



,u)=\og\ 



p- 1 A ft (u)dp+K 1 € 1 +K 2 € 2 \+C(log r) 2 






is a convex functional of log r in r <><//. where 



3 



tf=i(l+y2> #- 1 ; 



^UrlllogiR/roW+R^ro 1 . 



(2.2) 

(2.3) 
(2.4) 

(2.5) 



It was emphasized in [5] that the constants in [2.5] are usually not the optimal constants, which 
in any specific problem depend on the geometry of 2 and I) as well as the location of I\. 
An evaluation of the integral in (2.4) leads immediately to the inequality 



<r ^i + rolog (R/r Q )e 2 . 



I \ T p- l A p {u)dp-\r- 1 [ uhh 

I Jr Q " J w(r) 

Following the procedure indicated in [5] we are lead to the inequality 

r" 1 f u 2 ds+roh 1 + R^o%<2BM 2a - i) (K 1 ei+K 2 e 2 ) 6 . 

J co(r) 

The constant B is easily computed and in this case 

5=log (R/r)/log (B/r ). 
In the plate problem under consideration (1.5), (1.6) we now set 

Aw=u. 



(2.6) 



(2.7) 



(2.8) 



(2.9) 



Then (2.7) is valid with u replaced by Aw. Since we shall need it later we obtain now by 
integration the following inequality 

jj D {Aw)Hxdy+\ (f-rlWei+\ {r*-rl)Rh-£h 2 

<2B , r ,k„„,„ V / J -• (2.10) 



log [MmyriiK^+K^)] 



Note that for 



Kiei + K 2 e 2 <M 2 , 
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(2.11) 



which must hold if the estimates are to be useful, we obtain the cruder but simpler inequality 

ff(AtiOV8rfH-^ (2.12) 

We wish to find a bound for w 2 ds. But once such a bound is known then point-wise-- 

J u(r) 

bounds for w at points in D r follow at once; for we need merely introduce the biharmonic 
Green's function T(P,Q) which vanishes with its Laplacian on the circle of radius r and center 
at P , and has a fundamental singularity at a point P in D r . We then make use of the identity 

(2.13) 

where 2' denotes the portion of 2 which lies interior to K(r). An application of Schwarz's- 
inequality then yields the desired result. 

In order to obtain the bound for w 2 ds we introduce the following functional : 

J co(r) 

U(r,w)=\ogf [\-*\a p {w)+ {[^ (2.14) 



where 



m = w ^da, 712= WjiWjtdv, 

773= fl (Awydxdy, V4 =2BB 1 M^-^ ) (K l e + K 2 e 2 y 



(2.15) 



D Rl 

for 

r <r<7? 1 </?; £=log (i?//?i)/log (B/r ) (2.16) 

and 

P>1, £>P-1. (2.17) 

We demonstrate that for p and g satisfying (2.17) and the constants l h l 2 , h, h, and d suitably 
chosen, U is a convex function of r~ q . 

To this end we introduce the notation 

U (r, w) = log E(r, w) + dr- 2 * (2.18) 

(where E is the term in braces in (2.14)), and form 

q 2 E 2 j^- 2 =El (q+l-p)r 2 «+ 1 - p [A r (w)+ {{wAwdxdy]+r 2{ ^ l) - v ( f [|gradw| 2 

D r 

+wAw]ds\\-i^ +1 - p) \a t (w) + WwAwdxdyT+2q 2 C 1 E 2 . (2.19) 

D r 

As in [5] we obtain in a straightforward manner the inequality 

m 1 r 1l +m 2 r l2 <E~-^r- p ( w 2 ds-^(p-l) ff p-^^w'dxdy-h^-Urj^n.rj.+n^ (2.20) 

D r 
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whore 

-(j>-2) 






(2.21) 



2(/,+ l) l " 2 2(p+l) 

_ (p+2)r,-" _ r,-"-" 

^ /i+ ^+ir' n2 - h+ 2^+T) j 

The constants / x and l 2 are to be so chosen that rrii and m 2 are positive. 
By Green's theorem 

A r (w) + ( | wAwdxdy= ( to ^ <fe+ f w |^ &r. (2.22) 

Thus by applying Schwarz's inequality in (2.22) and using the fact that the arithmetic mean 
is not less than the geometric mean we obtain from (2.19) using (2.21), (2.7), and (2.9) 

q 2 E*j^ 2 >\ r~*f^ Ms | (g+i-^rit+i-ip [^(w)-| r- 2 jjw 2 dxdy-^- v ^\ 

+[^^jjr-«>+"w 2 dxdy+ m lVl + m*+^ + u] 

X-f (2+i-/»)/^ +| -' [A(a»)-^f- 2 JJw 2 cto/y-~ „,] 

+r 2«m>-/,rf |gradw| 2 cfe-^ f w 2 <fe-^-"j\ 

+2 2 2 6\ j^ / w2 ds+^^jjr-^+ l) wMxdy+l 3V3 +l iVi + m lVl + m^V- 

(2.23) 

where «i, 6 1; 6 2 , 63, and 6 4 are arbitrary positive constants. From the Rellich identity [6, 7], 
it follows that 

-L{^-^)}* rf J/'S'-**«X-'sSs^'5i»- 4 "i*. (W» 

Thus, for any 6x0, 

L{(S)*-®>^-<»>-^-V5»- (2.25) 

If b 2 U/(dr~ q ) 2 is to be non-negative, then the coefficient of the term A T {w) w 2 ds must be 

non-negative. After insertion of (2.25) in (2.23) we find that this coefficient is positive if 

g-7? + l-<5i>0. (2.26) 
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The coefficients of all other terms may be made positive by choosing d sufficiently large 
provided 

1>>1. (2.27) 

But (2.26) and (2.27) are just the imposed conditions (2.17). Hence for such a choice of the 
constants, U is a convex function of r~ Q . 
We now let 

8*=[(Ri/yo)- Q -m[(RiM- Q -l}. (2.28) 

It follows then that 

e w "^(r,w)<[6 c ^"^(Bi t w)] 1 - | *[«^o" ,, B(ro i w)r. (2.29) 

Using (2.20) we obtain the inequality 

\r~ v [ w 2 ds+&^ fw 2 r-^ l) dxdy+m l 7 ll +m 2 r j2 +l sm +hr J , 

<C i m^{^iVi+mri2Y*+G 2 Ml^ , {K 1 e 1 +K 2 e 2 y\ (2.30) 

where the constants C x and C 2 are determined as indicated in [5]. We have used (2.10) or 
(2.12) for 773 and the notation 

<5' = <5<5*. (2.31) 

By use of Schwarz's inequality and inequalities (2.7), (2.9), and (2.30), equation (2.13) then 
yields the desired pointwise bounds for w(P) at any point P in D r . 

In order to obtain the desired bounds in D— D r we first form from (2.13) w, iy Aw and Aw yi . 
For P a finite distance from the boundary of D r the differentiation may be carried out under 
the integral sign on the right and all of the resulting integrals converge. Thus using (2.7), 
(2.9), (2.13), (2.29) we obtain the pointwise bounds in D T not only for w(P), but also for the 
quantities |grad w\, Aw and grad Aw (in fact estimates for any derivative of w). If we now 
integrate these inequalities over a portion of a spherical surface oo(r } ), where r l <^r and each 
point on co(ri) is a finite distance from 2, then we may regard ca(r{) as a new Cauchy surface 
and proceed to define a new region D T1 . continuing as indicated in [5]. We are thus able to 
obtain pointwise bounds for w and its derivatives at any point on the interior of D. 

3. Analogous Viscous Flow Problems 

Instead of interpreting our biharmonic Cauchy problem as a plate problem we could also 
have considered v in (1) to be the stream function in the slow flow of an incompressible viscous 
fluid [4, p. 541] (here p = 0). However, in this problem the actual value of v is of no physical 
significance. The velocity components are expressed in terms of derivatives of v, and the 
physical quantities which we measure on S are the velocity components, the vorticity and the 
pressure gradient. In terms of v the vorticity is the Laplacian of v and the pressure gradient 
involves the derivatives of Av. 

Hence we are interested in estimating only the derivatives of v and not v itself. At the 
same time we would like for these bounds to involve only v\ 2) e 1} and e 2 , and not 771. In order to 

accomplish this we consider instead of (2.13) the identity 
w(P)-w(P) = - f [w-w(P)] £- (AT)ds- f Aw^fds 

~L {[«-«#>] I ^~fn AP+Aw d i~l <*"> r r*. t 3 - 1 ) 

where P is the point on S at which K(r ) is tangent to C. 
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For points away from the boundary of D T wc differentiate (3.1) carrying out the differenti- 
ation of T on the right. An application of Sehwarz's inequality yields the same types of terms 

as those encountered previously except for the terms [w— w(P) 2 ]ds and [w— w(P)] 2 da. 

Using the arguments employed in deriving (2.29) with w replaced by w—w(P), we may bound 
the integral over w(r) in terms of rj 2 , ei, e 2 , and [w— w(P)] 2 da. Let us now define X by 

r (¥)**> 

X=min Jo )r y (3.2) 

among continuously differentiate functions \j/ which vanish at the origin. The number X is 
the lowest eigenvalue for a vibrating string fixed at the origin and free at a=l. Thus 

X=tt 2 /4Z 2 . (3.3) 

We now let I be the shortest distance on 2 from P to either of the points of intersection of 
(a(r) and 2, and denote by L the length of 2. Thus 

l<L/2 (3.4) 

and since w—w(P) is clearly an admissible function \p in (2.32) — the origin being taken at 
P — we obtain using (3.3) and (3.4) 

f [w-w(P)]'da<L 2 ir- 2 I Q^) 2 da<Dir-\ 2 . (3.5) 

It is thus possible to obtain in D r the following bounds for the velocity components w tX using 

(2.7), (2.9), (2.30), (3.5), and (2.13): 

\w, x wJ 2 <( 1 [M\-^+C f 2 M 1 2 - 5, {K l e l + K 2 e 2 }^ (3.6) 

where C[ and C 2 are constants and all other terms are the same as those occurring on the right 
in (2.30). Note that as in [5] the constant M i? which is a bound for \w—w(P) | , may be replaced 
by a constant L u which involves a bound for the kinetic energy in D. 
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Selected Abstracts 

On Stokes flow about a torus, W. H. Pell and L. E. Payne, 

Mathematika 7, 78-92 {I960). 

This is the third of a series of papers on the Stokes (slow 
viscous) flow about axially symmetric bodies immersed in a 
uniform flow parallel to the axis of symmetry. The flow itself 
then has axial symmetry and is amenable to attack by the 
generalized axially symmetric potential theory of A. Wein- 
stein. In this particular paper the flow in multiple connected 
regions is considered and, in particular, the problem of flow 
about a torus is solved. The mathematical problem which 
characterizes the flow does not have a unique solution. 
Uniqueness is insured by the physical condition that the 
pressure be continuous throughout the flow field. 

Modern theories of materials, C. Truesdell, Trans. Soc. 
Rheology IV, 9-22 {I960). 

Theoretical work on rheological materials shows three distinct 
phases : 

1. Early work, to 1945. Simple models, mainly one-dimen- 
sional and linear, are constructed by analogy to discrete 
systems. Various coefficients are named. 

2. Intermediate period, 1945-1955. Ideal materials exhibit- 
ing nonlinear response are defined explicitly as exact mathe- 
matical realizations of a simple idea. Their properties are 
studied by means of general theorems and special exact 
solutions. 

3. Current phase, from 1955. The; most general constitutive 
equations consistent with the observed phenomena are 
determined. 

Work of the second and third kinds is summarized, with 
particular reference to the different roles of methods of 
invariance. 

Open problems are stated. 

Graphs of bivariate normal probabilities, M. Zelen and N. C. 
Severo, Ann. Math. Stat. 31, No. 3, 619-624 {Sept. 1960). 

Recently there has been much activity dealing with the tabu- 
lation of the bivariate normal probability integral. D. B. 
Owen [3], [4] has summarized many of the properties of the 
bivariate normal distribution function and tabulated an 
auxiliary function which enables one to calculate the bivariate 
normal probability integral. In addition, the National 
Bureau of Standards [1] has compiled extensive tables of the 
bivariate normal integral drawn from the works of K. Pearson, 
Evelyn Fix and J. Neyman, and H. H. Germond. In this 
same volume, D. B. Owen has contributed an extensive section 
on applications. 

It is the purpose of this paper to present three charts, which 
will enable one to easily compute the bivariate normal integral 
to a maximum error of 10~ 2 . This should be sufficient for 
most practical applications. Owen and Wiesen [5] have also 
presented charts with a similar objective; however, as pointed 
out below, we believe the charts presented here lend them- 
selves more easily to visual interpolation. Actually the 
motivation for these charts came from the Owen and Wiesen 
work. 

Comment on a paper of Mori on time-correlation expressions 
for transport properties, M. S. Green, Phys. Rev. 119, No. 3, 
829-830 {Aug. 1960). 

An auto-correlation expression given by Mori for the thermal 
conductivity of a fluid is shown to be only apparently different 
from an expression previously derived by the author. 



The Stokes flow about a spindle, W. II. Pell and L. E. Payne, 
Quart. Appl. Math. 18, 257-202 (I960). 

This is a continuation of earlier work by the same authors. 
The methods of the generalized axially symmetric potential 
theory of A. Weinstein and certain representation theorems 
of L. Payne for the solution of repeated operator equations 
are applied to the solution of the Stokes (slow) flow about a 
spindle-shaped body. The stream function of the flow is 
found and an expression is given for the drag of the body. 

Irrational power series, M. Newman, Prov. Am. Math. Soc. 
11, 699-802 {Oct. 1960). 

It is shown that if a. is a real number, g a non-constant poly- 
nomial, then 

2 g([na])x n 
ra=0 

is a rational function of x if and only if a is a rational number. 
The same statement is proved for the function 



X) x lnal . 



Topological derivation of the Mayer density series for the 
pressure of an imperfect gas, M. S. Green, J. Math. Phys. 1, 
No. 5, 391-394. (Sept -Oct. 1960). 

A new derivation of Mayer's classical density expansion for 
the pressure of an imperfect gas based on a classification of 
cluster graphs according to topological criteria is presented. 
The classification is a generalization of the classification of 
simple trees into trees with centers and trees with bicenters. 

Ensemble method in the theory of irreversibility, R. Zwanzig, 
J. Che in. Phys. 33, No. 5, 1338-1341 (Nov. 1960). 

We describe a new formulation of methods introduced in the 
theory of irreversibility by Van Hove and Prigogine, with the 
purpose of making their ideas easier to understand and to 
apply. The main tool in this reformulation is the use of 
projection operators in the Hilbert space of Gibbsian ensemble 
densities. Projection operators are used to separate an 
ensemble density into a ' 'relevant" part, needed for the 
calculation of mean values of specified observables, and the 
remaining "irrelevant" part. The relevant part is shown to 
satisfy a kinetic equation which is a generalization of Van 
Hove's "master equation to general order." Diagram 
summation methods are not used. The formalism is illus- 
trated by a new derivation of the Prigogine- Brout master 
equation for a classical weakly interacting system. 

On the theory of the critical point of a simple fluid, M. S. 

Green, ./. Chem. Phys. 33, No. 5, 1403-1409 (Nov. 1950). 

The consequences of a new system of integral equations for 
the theory of the critical point are discussed. Reasons are 
given for believing that the fundamental assumption of the 
Ornstein-Zernicke theory about the direct correlation function 
is incorrect. 
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